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PREPARATION OF THIS PAPER 

This document is based on and updates Fisheries Technical Paper No. S1 f Suppl, 1. 
It is one in a series of manuals and other documents prepared by the Fishery Resources and 
Environment Division to assist scientists at all levels ronoemed with marine research and 
its application to fisheries. It is concerned with one aspect of the use of tagging in 
marine research - that of the mathematical analysis of the results to provide estimates of 
Quantities (population numbers, mortalities, growth and dispersion) useful in further analysis 
of fish populations. It is planned to produce further reports dealing with other aspects of 
tagging - including the general purposes of tag/ring work, its advantage*! and diaad vantages - 
choice of tag; methods of capture and handling of fish; and methods of ensuring 0od reporting 
of recaptured fish. 
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Thin paper deals with the mathematical treatment of tagging data for estimating survival 
rates | population size, growth rates and movement parameters. It revises and updates an 
earlier paper with the same objective ( Jones, 1966), 

1 , Estimating Survival Rates 

One of the objects of tagging experiments is to estimate the survival rate of the indi- 
viduals in the tagged population. The final estimate may be expressed in terns of the pro- 
portion of individuals that survive* for one year (s) or as an instantaneous coefficient of 
total mortality (2), where 

S - exp (-Z). 

There are a number of methods for estimating survival rates and these may be conveni- 
ently divided into two groups: 

(a) Methods for estimating survival rate when tagging has been done on one occasion only. 

(b) Methods for estimating prurvival rate when tagging has been done on more than one 
occasion. 

1.1 TSert jjnat ing survival rate from a single tagging experiment 

If tagging is done on only one occasion, the survival rate has to be determined from 
the rate of decline of the tagged individuals. 

Recaptures may be made in one of two ways. One method is to take samples from the popu- 
lation on a number of occasions and to record the numbers bearing tags in each sample. This 
method is appropriate if it IB possible to tag so large a proportion of the population that 
the expectation of obtaining tag^od fish within subsequent samples is reasonably high. When 
dealing with many commercial fish populations, however, the proportion of the total stock 
that oan be tagged may be extremely small. As a result the expectation of obtaining tagged 
fish in subsequent samples is often also very small. In this situation, the experimenter 
may be dependent on commercial fishermen for the recapture of tagged individuals and, in 
this situation, the recaptures will be returned more or less continuously. 

Whichever method of recapture is used it is implicit in the methods used that the num- 
bers of recaptures can be expressed in units proportional to the total numbers of tagged 
individuals remaining in the population at successive times. There are various ways of doing 
this. To be precise, let 

n. be the number of individuals taken in a sample on the ith oocasion 
m. be the number of these found to bear tags 

y be some function of n^ that is assumed to be proportional to the number of tagged 
individuals in the population. 

There are at least three possible choices of the function y. . One is y. - m^/n^, which 

requires the assumptions that the proportion of tagged individuals in the sample is the same 
as in the total population, and that the size of the population does not change during the 
experiment. A second choice is y - n^A^t where f is the fishing effort expended in col- 
lecting the ith sample. If the sampling effort is the same on every occasion, the simpler 
form y. m. may be chosen* 

If the tagged individuals are recaptured continuously, then choices of y^ can be made 
In the same way. Thus, values of y may be calculated either from the ratio of numbers of 
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tagged individual* to total numbers of individuals caught over a period of time or they may 
be calculated from the number of tagged individuals taken per unit of fishing effort over 
some specified period of time, depending on the data that is available and the assumptions 
that are appropriate. 

There are a number of methods of estimation hued on the assumption that the survival 
rate is constant over the experimental period. Estimates of survival rate can then be 
obtained given successive values of y. taken either at equally spaced moments ir* time or over 

successive periods of equal duration. 

Table 1 shows four methods of estimation that can be used for estimating the survival 
rate from a single tagging experiment when the recaptures are made on discrete occasions. 
The values in the row m. refer to number of individuals recaptured on six occasions. In the 

row n. are given the total numbers of individuals obtained on each sampling occasion. The 

values in the row y have been calculated assuming that the proportion of tagged individuals 
in each sample is proportional to the number of tagged individuals ua the sea. For conveni- 
ence, these estimates have been multiplied by 10^. Thus, OB the second occasion 

93 - (140 x 10 4 )/15 000. 

Pour methods of determining the annual survival rate are given. The first method is 
due to Jackson (1939) and was referred to by him as his "positive method". His estimate for 
the survival rate is 

s " <y 2 + y 3* y k )/(y i +y 2 * W- 

Jackson also refers to an alternative formula proposed by W.M Stevens which is 

s - [(y 3 +y 4 + y k )/^i +* 2 + y *-2 

Chapman and Robson (1960) thought that these estimates suffered from the disadvantage that 
they could take values greater than unity and considered that a better estimate of the sur- 
vival rate was given by a modified estimate due to Heincke (1913) which is equivalent to 



Robson and Chapman (1961) consider methods of analysing catch curve** and their methods 
of estimating survival rate can be applied tc data such a thoae in Table 1. The procedure 
is as follows. 

First a total coded age (x) is calculated where 

X m (i - l)y for i - 1 f 2, ...J (the number of samples). 

Hie second step is to calculate the mean coded age (x) where 
x - X/t y t . 

The final step depends on the umber of tags still outstanding after the last (jth) sample 
has been collected. If, at the final sampling, it is supposed that the proportion of tagged 
individuals still at large is too large to be neglected then the estimate of mean coded age 
is equivalent to 
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J- 1 . J- 1 . 

x - r i s 1 / r s 1 . 

i - 1 i - 

Thio equation cannot be solved directly for S but estimates of S needed to satisfy any given 
value of x oan be determined for any value of J using Table 3 of Robaon and Chapman (1961). 
Thus, for the example in Table 1 f x - 1.5319 and_J - 6. Prom Robson and Chapman's Table 3 
it oan be seen that for 6 classes this value of x is satisfied by a value of r - 0.70. If 
sampling continues until the number of tagged individuals remaining in the sea is negligible 
then Robson and Chapman's method can be simplified. Suppose, for example, in Table 1 sampling 
had been continued on occasions 7 f 8 f 9 f etc., until no further tagged individuals were 
obtained. Suppose also that the entry given for occasion 6 represents the total numbers 
obtained over the 6th and all subsequent occasions. In that case an estimate of survival 
rate is given directly by 



+ X 



- y g ) 



Using the data in Table 1, for example, this would have given 
S - 553/(361 + 553 - 18) 
- 0.62. 

Table 2 illustrates the same methods of computation for the situation where values of 
y. have been obtained from the catch per unit effort of the tagged individuals on a number 
or ocoasiono. Each value of y, is given simply by the ratio of the number caught and the 
fishing effort (f.). Thus, for example, 

9.66 - 140/14.5. 

Table 3 chows similar calculations for the situation where individuals have been recap- 
tured continuouoly over a period of 4 years. The periods have been ooded from 1 to 4 &d the 
umbers recaptured hove been converted to numbers per unit effort to provide values of y,. 
Estimates using the methods of Jaokaon, Stevens, Heinoke, and Robson and Chapman oan then be 
obtained exactly as before. 

In addition to these methods, regression estimates oan also be made. These are baaed 
on the fact that, givon a constant mortality rate, the decline in the numbers of tagged indi- 
viduals can be oxpocted to approximate to an exponential curve. An exponential curve oan be 
made linear with reopoct to time by taking logarithms of the numbers alive at equal intervals. 
A regression of log cumbers against time should then provide an estimate of the instantaneous 
mortality coefficient. Z, directly. Tfcis method has been developed by Beverton and Holt 
(1956), Oulland (I969)aud Paulik (1963). 

Thus, suppose fishing effort is constant throughout the experiment and let y, be the 
number of tagged individuals recaptured during the ith period after tagging. Eacn period 
being of duration p (i.e., the ith period extends from tine (i - l)p to tine ip), it follows 
that 



1) p)) (1 - exp (- Z 

(i - 1, 2, 3 ..... ) 

where F is the instantaneous mortality rate due to fishing 
T is the number of tagged fish released 

Z it the instantaneous mortality rate due to fishing and all other causes, including 
tagging mortality subsequent to the tiae of tagging. 
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Thie equation can be made linear with respect to i by taking logarithm* of both sides* Thus 
log ( yi ) - -Z (i - 1)p + log (FT) + log ((1 - exp (-Z))/Z). 

This eh owe that by plotting the natural logarithms of y against the dummy variable i t a line 
of the slope -Zp is obtained (Beverton and Holt ? 1956; Jones 9 1956). Aw example ie shown 
in Table 4* Thin shows, for comparative purposae, the same numbers of individuals captured 
per unit effort as were used in Table 3* Natural logarithms of these numbers have been cal- 
culated and a regression of log (y ,) on i is found to give 



log (y i ) - 5-M7 - 0.430 i 

This equation shows that the estimate of Zp 0.430* Since the grouping period is one year 
p - 1 and hence Z - 0.430 also. 3 ie then obtained directly from the relationship 

S - exp (-Z) BO that 
S 0.6? 

Thi*; it can be noted, ifv eUaliar to values obtained ueing the method of Jackson, Stevens, 
azw? Robeorr and Chapman* 

Table 5 shows a similar r v.i;*> i <a using the saraa basic data but with the regression 
method of Gulland (1969)* Thie mthod differs from that of Beverton and Holt simply in that 
xhe coded periods are labelled 0,5, 1,% etc., instead of 1 f 2 f 3, etc, This method of 
coding is to simplify the estimation of F which will be discussed later* Again, the 
regression of log (y.) on i giveb a slope of -0*430. Hence Zp - 0,/f3f> and hence the estimate 
of S is 0.65 as was obtained in TV:. Vie 4* 



Table 6 shows the method de*M5r.bd L 4 PauliK (1963) baaed on the BJB* principle "but 
using a weighted regression technique. Pauiik ueae for weighte^ the values of y i /T f where T 
is the number initially tagged* In Table 6, T io t&kon as 1 QOO. Instead of values of 
log (y. ) the regression ie carried out uwing the logarithms of y./T. The numerical values 
are shown in the table and the weighted regression has a slope or -0.438. The estimate of S 
finally obtained is 0.645, eioilar to the value obtained by the previous methods. 

Some of the methods have been described with reference to the situation in which recap- 
tures have been made on discrete occasions ,, and other have been applied to situations where 
recaptures are made continuously* There seams no objection, in principle,, why both methods 
of analysis should not be applied to both situations. The distinctions above relate more to 
the way in which these methods have been developed by different authors than to restrictions 
in the way in which they can *b used. 

Pauiik (1963) compared the regression method shown In Table 6 with the i""rt"im likeli- 
hood method described by Robson and Chapman* He considers the regression method to be more 
useful for general application since it requires only easy computations regardless of the 
number of recovery periods involved. However, in the situation where the numbers of Indi- 
viduals returned in some of the time periods are snail or ssero, the regression method would 
not necessarily be satisfactory and in this situation the maximum likelihood estimates of 
Robson and Chapman are to be preferred. 

Another way of estimating a mortality rate from a single tagging experiment is described 
by Ottlland (1955). He considers the situation In which the times from liberation to recap- 
ture axe known for each fish recaptured and there are no tagged fish outstanding at the end 
of the experiment. He then assumes that the number of tagged fish declined exponentially 
with a constant mortality rate (z) and shows that T Is an unbiased estimate of 1/Z where t 
la the mean time of freedom of the recaptured individuals. To a first approximation there- 
fore 1/t is an estimate of Z. An example Is shown in Table 7* The basic data are the actual 
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days of freedom of the individual recaptures. In the example the tie have a mean of 63 day* 
which is equivalent to 0*17 years* The annual value of Z ia therefore equal to 5.9. 

Paulik (19&3) considered maximum likelihood estimates for the Oulland model when some 
tag* may still be outstanding ojad showed that an estimate of the antm^l survival rate (s) 
oan be obtained from the following equation: 



(S D ) - 



S/0 - 



where U is the duration of the recovery period. Provided t if? lees than half of the total 
recovery period this equation can be solved with the aid of a table given by Deemor and 
VotaM (1955, p. 501, Table 1). An example is given in Table 1. 

1.2 Estimating mortality from more than one tagging experiment 

A considerable number of theoretical papers have been written on the analysis of mul- 
tiple tagging and recapture experiments* By this is meant experiments in which Individuals 
are tagged on more than one occasion and sampled for recaptures on more than one occasion. 
In experiments in which there is a relatively high probability of recapture, multiple 
tagging leads to the possibility of recapturing and releasing the same individuals on more 
than one occasion. 

After the first sampling there may be individuals which have been recaptured on more 
than one occasion. If so, there are three ways of classifying the data: 

(a) A recaptured individual oan be classified with reference to all the occasions on 
which it has been previously tagged and released, 

(b) A recaptured individual oan be classified with reference only to the occasion on 
which it was first released. 

(o) A recaptured individual can be classified with reference only to the occasion on 
which it was last released* 

The first method has been considered by Jackson (1939, 1948) and Fisher and Ford (1947). 
Later workers have considered that this is a less efficient method than grouping individuals 
according to the occasion on which they were last released (Leslie and Chitty, 195 1 Moran 9 
1952). 

Figure 1 shows a general way of classifying the numbers returned. It is assumed that 
individuals have been tagged on sampling occasions 1 to 6. Recaptures have been made on 
sampling occasions 2 to 7. The numbers returned from each experiment and each sampling 
time are indicated by m. . (thus m p _ refers to individuals released on occasion 2 and recap- 

tured on occasion 5). If grouping method Bo. 2 is used this would refer to individuals 
first released on sampling occasion 2. If grouping method 3 is used however, this would 
refer to individuals last released on o ocas ion 2 (i.e., some of them could have been first 
released on occasion 1 but then recaptured and subsequently released on occasion 2). 

The numbers released on any occasion have also been classed in different ways by diffe- 
rent author*, toe nay is to count only the individuals that are newly tagged and released 
on each occasion. fRie other is to count the total number released on each occasion (i.e., 
imdividuals tagged for the first time plus any recaptures already bearing tags that are 
successfully released again). 

1.2.1 & Jol r~Sbr method 



The most general solution is one proposed independently ty both Jolly (19*5) and 
Seber (1965). The model is stochastic and, according to Cormack (1969) t the solutions 



Arrived at are so muoh simpler to apply than those obtained from the corresponding determi- 
nistic model that they should totally supersede them. 

In the Jolly-Seber model recaptures are grouped according to the last time an individual 
is released (method "o" above)* The numbers released are classified according to the total 
number released on eaoh occasion (i.e., newly tagged + previously tagged). The situation is 
illustrated diagrammatical ly in Figure 2. Tagging and capturing are both done during brief 
periods, between which are longer periods in which no recaptures are made and no tags are 
released* Consider, for example, the situation at the beginning of the third tagging period. 
At that moment there will be I ' tagged individuals in the sea. During the duration of the 
third tagging period nu of these will be recaptured and successfully released strain as part 
of the T, individuals released at the end of the third tagging period. This means that at 
the end of the third tagging period T. individuals will have been released and there will be 
(I. - O individuals in the sea bearing tags that were tagged prior to the beginning of the 
third period. It should be noted that I. will be smaller than I', because of mortality of 
fish in the sea during the course of the tagging period. Daring the remainder of their 
lives q~ out of these (l~ - nu) individuals will be recaptured. Similarly r^ out of the T-, 
individuals will be recaptured. It follows then that by simple proportion 



hence 1^ nu + q 3 T-jAy 

Similar considerations during the fourth experiment lead to the relationship 

'4 - "4 + 4 VV 

The proportion surviving from the end of the third experiment to the beginning of the fourth 
experiment is then given by 



S '3/4 ' X V (1 3 - -3 + 



3imilarly f the proportion surviving from the end of the third experiment to the end of the 
fourth experiment is given by 

S 3/4 ' V (I 3 * "3 * T 3> 
o that 

S 3/4 " r 3 (T 4 q 4 + r 4 *4 )/(r 4 T 3 (q 3 + '* 
In general it follows that the relationship for survival rate is given by 

t+i T t 

An example is given in Table 8 which shows the numbers tagged on four occasions and the 
numbers returned, classified according to the scheme in Figure 1 . It is assumed that all 
individuals captured that already bear tags are released again. Thus, for example 9 of the 
1 656 individuals tagged on occasion three f 156 of them were already bearing tags so that 
only 1 500 were freshly tagged. The first step in the calculation was to calculate the 
column totals (m.) and the row totals (r.). The next step was to calculate values of q noting 
that <h in all oases* Values of 8 were calculated using the formal a above. 
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1,2.2 Robson's method 

An alternative solution to the Jolly-Saber method ie one proposed by Robson (1963). In 
this method the recaptures and the numbers tagged are classed in a different way from that 
used in the Jolly-Seber model. They cure classed according to the occasion on which they were 
last released and the number tagged on each occasion refers only to the number newly tagged* 
It does not include any individual o that happen to be recaptured bearing tags and released 
again* The situation is illustrated diagrammatioally in Figure 3. 

Consider, for example y the situation at the end of the third tagging occasion* At that 
moment there will be I- tagged individuals in the sea. Of these q^ will be captured subse- 
quently* At the same time, T^ individuals will be tagged and of these r. will be captured 
subsequent ly. It follows by simple proportion that 



Similarly for the fourth experiment it follows that 



proportion surviving from the end of the third to the end of the fourth experiment is 
then given by 

S 3/4 - V (I 3 + T 3 ) 



" T 4 q 4 V (T 3 r 4 (q 3 +T l^' 
A gmeral formula than follows dirotly fron thiit 

S t/t+1 - T t+1 t-H r A r t+1 ( <r + p t ) ' 

An example of the method is given in Table 9* It is assumed that the recaptures have 
bean taken during discrete intervals at the tagging times* Hote that q, - and that other 
values of q are calculated as indicated in Figure 3. The estimates of the survival rate have 
been made using the formula above* 

It should be noted that Robson's formula also lends itself to the situation in which the 
recaptures are made continuously and independently of the times when tagged fish happen to 
be released* 

Robson's formula is equivalent to that of Jolly and Seber if the terms n^ happen to 
equal eero. Tfcis would be so if the probability of recapturing a taggod individual during 
the course of tagging happened to be zero, as is liable to be the oaae when dealing with 
large populations* 

An example is given in Table 10* It is assumed that the recaptures are made continu- 
ously, and that these have been grouped into periods corresponding with the periods in between 
each tagging experiment* 

1*2*3 The triple patch method 

Olhis is a apeoial oaae, dealing with the multiple tagging and recapture situation when 
there is sampling on three occasions only. Cta the first two occasions individuals are tagged; 
on the second aad third occasions recaptures are recorded. There are two solutions depending 
em the method of grouping adopted* 
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Bailey (195^) grouped the recaptures according to the tiine of first marking. He 
obtained maxiuum likelihood estimates of survival rate from the end of the first sampling to 
the end of the second sampling occasion. Hie equation is 

S 1/2 ' T 2 

In this situation BU-> is equivalent to Robson's q ? . Also m ^ is equivalent to Bobson's r~* 

The result is equivalent therefore to Robeon'e result when i 1 (note that q , m o). An 
example is given In Table 11, baaed on data taken from Table 9. The same va3ne of 0.60 is 
obtained for S / using Bailey's equation as was obtained using Robson's general equation. 

The triple oat oh method waa also considered by Leslie and Ghitty (1951 ). They grouped 
recaptures according tc the interval when they were last released and the estimate for sur- 
vival rate from the end of the first sampling occasion to the beginning of the second is 
given bry 

S 1/2" ( 

This result is equivalent to JoHy and Sebor'e in the situation where i 1 f provided all 

of the individuals recaptured ora the socond occasion are re-trained unharmed to the population* 

An example is given in Table 12 uoln# data from Table 8* The saai* value of 0.59 for S 1 / 

is obtained using Leslie artf Chitty's equation as was obtained by Jolly and Seber's general 
equation. 



1.3 The estimation of fishin/^inorl^ 



Tagging experiments not only rtfeble estimates to be made of annual mortality rates but, 
In principle at least, also enable tbeB^ awtirateg to be split into component* due to fishing 
mortality and due to other 



In the situation wherr* only one ta^cg^ug haa been done, the regression methods of Beverton 
and Holt, Oulland and Poulik all provide estimates of P t the instantaneous rate of mortality 
due to fishing* I-a the oese of the Beverton ciad Holt ro&ression f the estimate of P is given 

tor 

P - Z(exp (I o - Z P ))/(T(1 - r 

where I io the intercept. 
Por the GullaDd method the esti^ite 
P - (e*p ( 



For the method decrilftd by Pmilik this ^tinata is 
P - 2 (exp (I o - 7, P ))/(1 - e-i> (-Zp)). 

Examples are given in Tables 13-15 using data from Tables 4-7. Por all the examples a 
value is required for the initial cumber tagged and this has been assumed to be 1 OCX) indi- 
viduals in all cases. Table 13 shows the application of the Beverton and Holt, azvl Paulik 
methods to the data ffom Tables 4 and 6* Prom Table 4 the intercept I - 5347 The value 
0.17 of P then follows from the formula. Prcii Table 6 the intercept iSt -1.539; substitution 
of this in the appropriate equation also gives a value of P 0.17. (it should be noted 
that the difference between TWO intercepts is that in the weighted regression the number* 
recaptured were all divided by 1 000. Apart from statistical variations in the fitted 
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regressions, the differences between the two intercepts is therefoie equivalent to the 
natural logarithm of 1 000 which is 6*99*) 

An estimate of mortality rate can also be made in the cane of those experiments in which 
the exact liberation times for each individual ar known. Th* Oilland (1955) formula for 
estimating F is 



F - 

where m is the total number recaptured, 



Chapman (1961) showed that for this wodijl no minimum variance unbiased estimate of F exists 
but that the following estimate of F is almost unbiased: 

p. (,-1 )/(). 

An example is shown in Table 15 ueing data from Table 7 azad assuming that the number of fiah 
tagged ir the first instance* was 1 000 individual n. In thin example the values of F obtained 
using the Oulland and the Chapmnai formula* acre both equal to 0,12. 

An estimate of X which is equivalent to the mortality rate due to all oauetas other than 
fishing is equal to 7, - F. In the example in Table 15 thie is *qual to 5.9 - 0.12 - 5.78. 

When dealing with multiple tagging experiments there are two alternatives to consider. 
Tables 8 and 9 deal with the situation in which sampling for recaptures is done on the same 
occasions in which fioh are captured for tagging. In this situation vhe rat* of capture on 
the ith occasion (p.) is given "by 



To a first approximation I' - I. provided the duration of each taking experiment is negli- 
gible in illation to the total lime in which tagged fish or a tagged end recaptured. An esti- 
mate of p. is then given by ^./I.. Examples are given in Tables 16A and B using data from 
the examples in Tables 8 aBd y. jEstimatee of I. have boen taken from the basic formulae 
developed ia Figures 2 a&d 3* 

Estimates of fishing mortality F (as distinguished from the rate of capture p. above) 
can also b made in the case of maltiple tagging experiments. Table 1? f for example, illus- 
trates a way of doing this unin# the data front Table 10. In this example it has been assumed 
that the recaptures have b<sn mr.de oont inuouftly durii# the periods between tagging experiment*. 
Thtu) 9 for any one period it is possible to write 

- F N (1-S)/S 

where m is the number caught and H is the number of tagged individuals in the sea at the 
beginning of a recapture period (and hence at the end of a tagging period). Estimates of H 
at the end of aay recapture period are equivalent to T^ + I i (see Figure 3) and this i* 
given by 



Given estimates of I., n. and Z., it i* a simple matter to determine the value* of F. 
Detail* are siren liTTabt* 17* *aaio data have been taken frca Table 10. Hie valve* of 1 
hare ben calculated using the formula above. Value* of ^ have been calculated frow the 
relationship 

Z - -log (S). 
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Estimates of P have then been obtained from the relationship 
P - Z 



1.4 Some factors that can influence the estimates of survival rate obtained from tagging 
experiments 

There are a number of systematic errors that can influence tagging estimates and these 
hccve been classes as types A, B and C by Ricker (195*0 T^i* underlying assumption of all 
methods of analysis is that in the constant parameter situation the numbers of tagged fish 
will tend to decline exponentially with time. This means that the logarithm of the numbers 
alive at any moment should be linearly related to time as shown by the continuous lines in 
the diagrams in Figure 4. l*y type A errors are meant those errors that tend to influence 
th intercept but not the slope of the regression. This, for example, can happen if there 
is an initial mortality of tagged fish due to the act of tagging itself so that the number 
of live individuals liberated is smaller than the nunfber actually tagged in the first 
instance, lype A errors can also occur if there ie incomplete reporting of tags. In parti- 
cular, if a constant proportion of the individuals recaptured are reported the effect will 
be to move the regression linos down, influencing the intercept but not necessarily the 
slope. If the elope of the regreoeion ie not effected, estimates of survival rate and total 
instantaneous mortality rate (Z) are not affected. The effect of an underestimate of the 
intercept, however, is to underestimate estimates of fishing mortality rate (P). 

ErroirB of types B and C tend to influence the slope of the regression but not neces- 
sarily the intercept. The situation is illustrated diagrammatic ally in Figure 4. l^rpe B 
errors occur if the act of marking causes marked fish to die at a greater rate than 
unmarked fish throu^iout the entire course of the experiment, Tliey could also occur if there 
was pz-ogressive loss of interest in the experiment by fishermen so that the proportion of 
tags returned tended to decline during the experiment * Type C errors are those that cause 
the Blope of the regression to be too great due to an emigration of marked fish from the 
experimental area* 

The effect of errors of types B and C is to underlet imate the survival rate and over- 
estimate the total instantaneous mortality rate (Z). Errors of types B and C do not neces- 
sarily affect estimates of fiching mortality rate, however, providing they are made from the 
formulae of Qulland and Chapman which are not dependent on the estimates of Z. 

2. The Determination of Population Size from Marking Experiment o 

Marking experiments can be used for determining the site of fish populations* In prin- 
ciple the methods exe esrtremoly simple and are based on some variation on a method of esti- 
mation known as the Pet or sen method. 

2.1 The Pet era en method 

The Peterson method consists of tagging individuals on one occasion and sampling for 
recaptures on a single subsequent occasion* Suppose, for exanple, that there is a population 
of H fish and that T of these are tagged* The simplest situation is when the tagged fish are 
distributed at random throughout the population, or when the sampling effort is distributed 
at random. If th'a proportion of tagged fish in a sample should be the same as the proportion 
of tagged fish in the population, i.e., if there are m tagged fish in a sample of sime n, 
then 

m T 

is expected to equal o, 

00 that on rearranging terms one gets 
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In the literature, the estimation formula does not always appear in this simple form 
because it is biased and needs to be corrected. The cause of the bias ie toat w f the number 
of marked fish in a sample of size n f is a random variable. Over many samples, the observed 
values of m/a will indeed average T/N, but it does not follow that the observed values of 
the reciprocal n/m f which is what appears in the estimation formula, will average H/T. 
Instead, their average will exceed R/T because, in technical terms, the expectation of the 
reciprocal always exceeds the reciprocal of the expectation. The simple formula, therefore, 
consistently overestimates population size, 

A number of modifications of the basic equation are fflven in the literatv^e to correct 
this bias* These are summarised in Table 18* Eouations (A) - (c) are used when the Campling 
is "direct 11 ; i.e., sampling until a sample size (n) has been obtained* (The eamplo sice may 
be predetermined but there is no theoretical requirement to fix n before sampling starts, 
and in practice it is usually impossible to control sample size anyway.) An alternative pro- 
cedure, known as "inverse" sampling, is to continue sampling until a predetermined number of 
tagged animals (m) has been recovered. The appropriate equation for this situation ie given 
in (D) of Table 18* 



Examples of population estimation usijog the Petereen method ar* given by Rioker (1958), 
Hancock (1963), Le Cren and Kipling (1963), Sato (1938) and Simpson (1963). 

The reason for the extreme simplicity of the Pet or pen method lies in the assumptions on 
which it is baaed. The fact that either the fish or the sampling effort should be randomly 
distributed has already been referred to. Another important assumption is that the tagged 
and untagged fish are equally liable to capture and, if thies does not hold, population esti- 
mates are likely to be biased (Carlander and Lewie, 1948; Junge, 1963). Anderson and Bagge 
(1963) describe how plaice tagged with Peteroen discs were caught up in the material of 
nesheis through which they would otherwise have escaped. Their chances of capture were thereby 
increased; to estimate the increase further experiment e had to be conducted using a diffe- 
rent type of tag. 

2.1.1 TO number of fish that have to be tagged 

The numbers of fish that have to be tagged and the siae of the sample that has to be 
subsequently taken are directly related to the siee of the population. 

According to Robson and Regier (1964)1 to he sure of getting an unbiased estimate, it 
is necessary to take a sample large enough to be certain of recapturing at least one tagged 
fish} i.e., if T fish are ta#ff&d out of a population of H flab, then the subsequent sample 
(n) should be not less than -T, They show that if this oondition Is not met the resultant 
estimate of H will be biaced by approximately 100 exp (-(T + 1)(n + 1)/K) percent. In order 
to ensure a negligible bias. Robson and Regier recommend that the product of T x n should 
exceed the population sine (if) \iy a factor of at least 3 or 4* Robson and Regier give 
graphs for estimating how many fish should be tagged and how many fish should be sampled 
subsequently in order to estimate population sine with various degrees of precision. 

According to Cormaok (1969) the distribution of nT/ is far from symmetrioal , so that 
confidence limits oanaot be set to the estimate of population size in the normal way. As a 
way out of this difficulty he suggests that confidence limits should be set to m/nT, since 
the distribution of this is more likely to be symmetrical . The reciprocals of these limits 
should then provide confidence limits that are lees biased. 

Davis (1964) gives graphs for setting confidence limits directly to population sice for 
different values cf population sime, number tagged and number subsequently sampled. 

2.2 Other methods 

Im practice it My not be possible to take a sub-sa*ple of the population just after 
tagging. The longer the interval that elapses, however, the smaller is the likelihood that 
Petersen'0 equation in its simple form will still be applicable. 
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T 
The assumption that c is constant throughout the period in which recaptures are Made 

is implicit in the Pet ere en equation. Ideally therefore, the Petersen method should "be used 
before any mortality of tagged or untagged fish has had time to occur. The method is not 
necessarily restricted in this way, however, provided the mortality and/or emigration rate 
of tagged and un tagged fish is the same. For example, suppose that after a certain period 
only a proportion, p, of the original population of H individuals remains alive within the 
sampling area* The losses may be due to mortality or emigration from the sampling area or 
to both. The population size is then Hp. Provided the proportion of tagged rish that 

remains is also p, it follows that the ratio of tagged to un tagged fish will be 7^ which is 

T P 

still equivalent to =. If this is so f the basic equation oould still be valid for estimating 

population size, even though tagged fish may have been recaptured over a considerable period 
of time, e.g., when tagged fish are recaptured by commercial fishermen over a period of 
several years. 

If the recapture of tagged fish occurs over a long enough period, some account must be 
taken of changes in population size due to recruitment or immigration. If this happens, H, 

T 
but not T, may actually increase so that the ratio j: will then no longer be constant. 

Instead it will decrease. The simplest situation is one in which the recruits, or immigrants, 
can be distinguished by their sise or some other feature, from the original population. If 
that is BO T they may be excluded from the sample of n individuals and then the basic Petersen 
method can still be applied. 

If recruits or immigrants cannot be excluded, an alternative procedure due to Parker 
0955) rosy he adopted. This method is based on the assumption that a decline in the pro- 



portion of tagged fish in the population will be reflected in a decline in the proportion of 
tagged fish in the samples; i.e., if jj declines, so also will , where T is the number of 

tagged individuals in the population subsequent to tagging. Parker v s method then consists 
of plotting the ratio against time. A line drawn through these points will have a nega- 
tive slope, and if it is extrapolated back to zero time it will have an intercept that rep- 

* m 

resents the value of at the time of tagging. This, of course, will represent before any 

n n 

recruitment or immigration has had time to occur and so will be the value of use in the 
Petersen equation. This method has been used by Hancock (1963) for estimating the size of 
whelk populations. 

In a later paper Parker (1963) describes a method of estimating the rate of recruitment 
by using a model in which the rate of recruitment is assumed constant throughout the experi- 
ment. A system of equations is derived for estimating the population size, the rate of 
mortality and the rate of recruitment. These cannot be solved algebraically but a solution 
can be obtained by using the appropriate numerical technique. 

Jolly and Seber give an estimate for population size which is identical to the Petersen 
equation; i.e. 

H i - n i Vv 

The importance of the Jolly-Seber method, however, is that an aquation is also derived for 
estimating the number of tagged individuals in the sea at the time of the ith sampling; i.e., 
they show that 



so that by substitution 
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"l " (n i "i T l /m i r i ) + V 

An example is given in Table 19, using data from Table 16A, as well as a set of values of n. 
provided for the purpose of the example. Since values of I, have already been calculated, 
the estimates of population siae follow directly from the Peterson equation in its simplest 
form* Vote that I. is here assumed equal to I'. Similar estimates oan also be obtained 
using values of If obtained by Robson's method (Table 161.)* Estimates of H then follow 
directly from the equation above. 

Sohnabel (1938) considered an extension of the basic Pe1.*rnen equation 1r thf special 
situation where there is no mortality or recruitment to the population. The method is a 
simple extension of the Petersen method to a series of samplings for tagging and recapture* 
Sohnabel gives the following maximum likelihood equation for estimating the population sice H. 

(n, - 0^)1 



v -h 



It is assumed that n./N is sufficiently small to ignore the complications of sampling without 
replacement (say, leis than 0.1 According to Seber (1973)). 



equation can be solved for N, the population size by numerical methods (Seber f 1973; 
De Lury, 195^ ) An algebraic solution is possible too if various assumptions are made. 
Sohnabel considered in particular the oaae where the number of tagged fish in the sea at the 
time of the ith sampling happens to be negligible compared with the population size N. If 
this is so t the equation can be simplified to give the following estimate of NX 




Chapman (1952) suggested that a better estimate would be 




whilst Schumacher and Eschmeyer (1943) suggested that more weight chould be given to the 
numbers reoaptmred from large markings than from email by writing 




It has already been noted that these methods are based on the assumption of sero recruitment 
and mortality, and these are conditions that may only rarely be valid in practice. 

In the case of the "Triple Catch" method, Bailey also derived estimates of p, the sice 
of the population on the second tagging occasion, and R^, the proportionate increase in the 
population between the second and third occasions due to recruitment and/or immigration. 
These estimates are: 

*2 - T 2 *2 "l3 /m l2 "23 
H 2 - *3 "12/ B 2 "13 
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3 Rstimating Growth from Marking Experiments 

Marking re milts can provide valuable data for measuring growth y providing the marking 
process does not interfere with growth. This can be especially important for those species 
that cannot be aged or that oan only be aged with difficulty. 

3.1 Pitting a growth curve to data representing equal time intervals 

The basic requirements are, first to construct an age/length or age/weight curve, and 
second to determine the parameters of a curve that fits the data well. For the second pur- 
pose, the von Bertalanffy curve is now widely used (von Bertalanffy, 1938, 1949) and details 
and examples of methods of fitting this curve are given by Qulland (1969). 

32 Pitting a growth curve to data representing unequal time intervals 

The methods described by Qulland (1969) are suited to the estimation of von Bertalanffy 1 s 
growth parameters from growth data collected over equal time intervals* The growth incre- 
ments of marked fish, however, extend over varying lengths of time and some modification of 
the usual method is required to deal with these data. A way of doing this has been des- 
cribed by Oulland and Holt (1959). 

The basic von Bertalanffy equation is 



l t - LJ1 - .- - *o>) 



where 1 - the length at age t and L^ , K and t are the growth parameters. After some time 
interval a, the age will be t 4- a and the length will be 

+ a - t K 



The growth increment can then be expressed by 



Since the increments of marked fish will have occurred over varying time intervals they 
oan be standardized to a first approximation by expressing them in terms of increments per 
unit time. What is required is a quantity 

y <*t + a - V/a 
This is equal to 



The next step is to relate this quantity to the average length during the growing period, 
i.e., to the length \ + a * \. It this is called x, 



Rearranging terms gives 
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- * ) 



Next, substituting the expression for y above f gives 



Finally, putting *Ca - b gives 
tanh b 



1 
where tanh b - ( ~ 



e~ 2b 



or, alternatively 



From this equation it is clear that if y ^^ ^ oan ba plotted against x, the result 
be a straight line with a slope of -K and an intercept of K L^. 

An example is given in Table 20. The data in the first three oolxnnns are the recapture 
details for five fish. The values of y f the increment per unit time, and x, the mid point 
of each growing period, are calculated as shown* The next step is to plot y against x. The 
points will be found to li on a straight line with a slope of -0.2. The first estimate of 
K which oan be celled K^ is therefore 0.2. Using this value for K, b 1 can be calculated for 

each fish and the values of . i^ ; oan be calculated or looked up in the table given by 

1 y b i 

Oulland and Holt (1959)* In the last column of Table 19 are shown the values of 7 r~~ 
for each fish. tanh b 1 

b 1 

In this example, the values of . , , are so close to vnity that the values in the 

~tann D_. 

last column of Table 19 are effectively the same as those of y. If they had not been tht> 

y b-i 

same, the procedure would have been to plot , . > against x to obtain a new line with a 

tann D- 

new slope. The slope of this line would have given a second estimate of K(lO and the whol 
process oould then have been repeated until the estimates of X remained unchanged. In 

practice, it should be noted that for values of b up to about 0.4r the values of ^^ b 

remain close to unity so that a value of K can usually be obtained, as in this example, 
directly from the first plot of y against x. 

When y is plotted against x using the above data, the line is seen to out the x axis at 
a length of 90 cm* In other words, at a length of 90 cm the growth increment is ftero, so 
that 90 cm is the required estimate of L^ . 

The parameter t cannot be estimated from the data given in Table 19 alone. It is also 
neoessary to know the age of the fish at each of the 10 lengths recorded. If this is known, 
10 values of t may be determined, one from each length, aad the values averaged to give a 
single value of t . Details are firm *y Qollaad (1969). 
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4* The Determination of Movement Parameters from Marking E. er iment B 

The simplest way of demonstrating the pattern of movement of a group of recaptured 
marked fish is piotorially. The position of recapture of eaon fish can be plotted on a chart 
in relation to its position of liberation and its mean direction and distance of travel can 
be indicated* The number of days at liberty can be indicated too ag&i2urt each recapture 
position if desirod. 

To do this for each fish individually it* quite straight forward but there remains the 
problem of how to combine the values from individual fish so as to estimate the mean direction 
and velocity of travel 343 well as a coefficient of dispersion* The values required are? 

Y the mean direction of movement 

V the mean velocity in this direction 

2 
a the mean square dispersion coefficient (Skellam 1931)* 

41 Direction of movement (?) 

A way cf determining the mean direction of movement can best be appreciated from 
Figure 5 Consider the case of a single fish that has been liberated at the centre of some 
conveniently chooen sot of coordinates (A) and thai its subsequently recaptured at position 
B. The position of recapture can be represent txi within the coordinate system by giving it 
a displacement r and a direction 0. This method of representation is convenient in practice 
since r and are both valueo that crua be meaningfully assigned to any recapture position* 

For any one fish thin method of representation is excellent but for more than one it 
is useful to transform to cartesian coordinates. If the axes of the reference system are 
labelled x and y it is easily owon that the displacement of a recaptured fish from the origin 
becomes 

r sin 8 along the x aid.*, and 
r cos Q along the y axis 
That this is so can be demonstrated from the right angled triaagle (ABC) in Figure 6. Sinoe 

f - sin S 

it follows that T - r sin O. Similarly, 
* - cos Q 

and therefore y - r cos 0. If there is more than one fish tbeir total displacement along 
the x axis becomes 

r^ sin Q^ + r^ sin + - Er^ sin Q A 

Similarly, the total displacement along the y axis will be giTen by rr. oos 0. 

i * x 

These two displacement s f r sin d and rr cos can be used to locate what is effectively 
the net displacement of the group as a whole. This is done by choosing a point B such that 
th values of x and y (the coordinates of B) are givan by x - jjr sin and y r ccs (see 
Figure 7), Then, let R be the overall group displacement IB and let its component of 
direction be given by the angle j (Figure 7). 
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It follows that since 

- - tan T 



An example is given in Table 21 , wherein 
r sin Q - 1.23 
r oos C - -26.33 

tan Y - 

A tangent of 4O.0467 indicates an angle of 2.?*. Therefore 
T . 180 ~ 2.7 - 177.3. 

(Bote: Since Er sin is positive and Er cos is negative it follows that * lies in the 
quadrant from 90 to 180.) 

4.2 Mean velocity (?) 

In Figure 7 the overall displacement of the group ae a whole is denoted b^ R. Kbw t the 
total time that will have elapsed will be the sum of the periods of liberty of the individual 
fish. If these are t 1 , t ? , etc., the total time of liberty can be denoted by t. It follow! 
then that an estimate of 7, the displacement velocity of the group as a whole is given by 

S ; that is, by 



-/ f & sin O) 2 + ( 2r oos Q) 2 
* rt 

Determine a value of V from the data in Table 21 




4.3 M*an s 



is one more coefficient that is needed in order to describe the distribution of 
recaptures more fully. Tfcat is, a parameter is needed to measure, in some way, the extent 
to which individual fish move independently of one another. A high value of V and a low 
value of a 2 , for instance, would signify a real directional movement with a low rats of dis- 
persion about the group mean. A low value of V and a high value of a 2 would be suggestive 
ef random movement with comparatively little directional movement. 

As a basis for the estimation of a suitable dispersion ooeffioient the theoretical 
treatment >y Stalls* (1951 ) * Bavarton and Bolt (195?) of the physical theory of heat 
ocajdnoticn has been employed. Details of the theory have already been published for the 
situation where the direction of movement (T) is known beforehand (Jones, 1959)* I* is only 
necessary hare to rafer to the result obtained there and than to consider ita axtanaion to 
the oaaa where T la unknown. 
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In the case where Y is known 
( r cos Q' 



2 IJYr 2 
n (^ t 



In this formula it is assumed that each value of is measured from the actual direction of 
movement. Since in practice, it is convenient to meamire from some direction such as true 
north, it is (0 - Y) rather than that should be used in the equation above; i.e., if 
refers to the displacement from true north then (0 Y) refers to the displacement frorr* the 
actual direction of movement. Thus, 



2 1 r 2 rrr COB (0 - T )f 
n V- T ~ t 



l 
j 



Since Y niay be unknown it is desirable to eliminate Y from the equation. If this is 
done, (Jones, 1966) the equation for a^ simplifies to: 

o o o 

( r cos Q)*" + ( r sin Q) 



2 
Example; Determine a value of a from the data in Table 21 



- 18 km 2 /day 

2 
4.31 Significance of a 

The mean square coefficient of dispersion (a ) is really a compound parameter resulting 
from the product of c , the nean distance travelled by each fish between each change of 
direction and v, the average swimming speed of each fish. That is f 

2 e e 2 

a e.v or e. ~ 

where w is the mean time required to swim a distance e* The fact that the numerator is the 
square of a distance explains why a 2 is expressed in ^inito such a kn^/day. T^ie fact that 
it is a compound parameter results from th fact that a given rate of dispersion could arise 
because the fish swim elowly with infrequent changes of direction or rapidly with frequent 
changes. Without further information it is not possible to discriminate between these alter- 
natives. 

A numerical value of a 2 such as 4 Icn 2 /day caa be employed in the following ways. First 
if 4.0 a* > c.v f a table of * and v ca ba made up to satisfy this relationship. This is 
done in Table 22. If on uppar limit oan be set to the likoly nean wimraing speed v f then a 
lower limit is automatically set to the value of c. 

Second, a numerical value of a 2 can be substituted in the following equation to determine 
the proportion of fish expected to lie within a given distance of the marking position after 
a given time* The equation is 

/ 2 
P (R, t) - 1 - erp I- ~- 
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whara P (R t t) aignifiaa tha proportion of fiah lying within a cirola of radiua R f aftar 
time t. 

Nora generally, the equation 

, R 1 2 \ / R 2 2 

P (IL, R , t) - .IP - -i- - xp I- -f- 

v V a % 

can be need to determine the proportion of fiah lying between two oirolea with radii R 1 and 
H reapaotively. 
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Estimation of survival rate from a single tagging expeilraent 
with recaptures on discrete occasions - I 



i 


1 


2 


3 


4 


5 


6 


. 


120 


140 


16 


35 


40 


?9 


n i 


10 OOO 


'< r , OOO 


\ 000 


fl (XX) 


12 000 


16 OOO 



120 



44 



13 



Jackson 

S - (93 4- 53 * 4^ 4- 33 4- l8)/( 120 4- 93 + 53 4- 44 4- 33) 

- 0.70 

Stevens i 

S - [(53 4- 44 4 33 4- 18)/(120 4 93 4- 53 4- 44) T 

- 0.69 

Helnok* 

S - (93 4- 53 -f 44 * 3j H- 18)/(120 -f 93 4- 53 4- 44 4- 33 -^ 18) 

- 0*67 

Rob son and Chapman 

E t (i - 0^ - 93(1) + 53(2) 4- 44(3) + 33(4) 4- 18(5) - 553 
iy i 361 

i - 553/361 ~ 1.-5319 
S * 0.70 
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Table 2 

Estimation of survival rate from a single tagging experiment 
with recaptures on discrete occasions II 



i 


1 2 


345 


6 


a. 120 140 


16 35 40 


29 


f i 


10 14.5 


3.2 7.9 12.3 


15.8 


y i 12 


.00 9.66 


5-00 4.43 3.25 


1.84 


Jackson 
S - 




(9.66 + 5.00 + 
(12.00 4- 9.66 + 
0.70 


4.43 4- 3.25 4- 1.84)/ 
5.00 + 4.43 4- 3.25) 




Stevens 

S m 


[(5.00 4- 4.43 4- 3.25 4- 1.84)/(12.00 4- 9.66 4- 5.00 
0.68 


+ 4-43)]*" 


Heinoke 
S - 




(9.66 4- 5.00 4- 
(12.00 4- 9.66 4 
0.67 


4.43 + 3.25 4- 1.84)/ 

5.00 4- 4*43 4- 3.25 4- 1.84) 




Rob son and Chapman 
Ei (i_l) yi 9.66(1) 4 


5.00(2) + 4.43(3) + 3.25(4) + 1.& 


4(5) - 55-15 


y t 


- 36.18 






X 


- 55.15/36.18 








1.5243 






S 


- 0.70 







FIBS T153 



Table 3 

Estimation of survival rate from a single tugging experiment 
with continuous recapture* ~ I 



Months after release 
Coded period 



1-12 

1 

208 

1.5 

139 



13-2/1 
2 

182 
2.0 

91 



2>-36 

3 

156 

3.0 

52 



Jack*) on 

S - (91 + 52 + 40)/(139 + 91 + 52) - 0.65 

Stevens * 

S - [(52 + 40)/( 139 + 91)]* 0.63 

Be in eke 

S . (91 + 52 + 40)/( 139 + 91 + 52 + 40) . 0.57 

Robson and Chapman 

2(1-1) y - 91 + 52(2) + 40(3) - 315 
J5r i - 322 

x . 315/322 - 0.9783 

S - 0.64 



37-4^ 
4 

32 
0.8 

40 



Table 4 

Estimation of survival rate from a single tagging experiment 
with continuous recaptures II 

Regression method of Beverton and Bolt 



Months after release 
Coded period (i) 

number recaptured per 
unit effort (y. ) 



1-12 


13-24 


25-36 


37-48 


1 


2 


3 


4 


139 
4.935 


91 
4-512 


52 

3.952 


40 
3.690 



Regression of log (y. ) on i gives 

log (y t ) 5-347 - 0.430i 
Benoe Zp m 0.430 
Hence Z - 0*430 

S - ezp - Z - 0.65 
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Table * 

Estimation of survival rate from a single tagging experiment 
with continuous recaptures - III 

Regression method of Oulland ( 1969) 



Months after release 1-12 13-24 25-36 


37-48 


Coded period (i) 0.5 1.5 2.5 


3.5 


Dumber recaptured per 
unit effort (y ) 139 91 52 


40 


Log (y t ) 4.935 4.512 3*952 


3.690 


Regression of log (y.) gives 
log (y i ) - 5-132 - 0.430i 




Hence Zp - 0.430 




Hence Z - 0.430 




S exp - Z - 0.65 





Table 6 

Estimation of survival rate from a single tagging experiment 
with continuous recaptures IV 

Weighted regression method of Paulik ( 1963) 



Months after release 
Coded period (i) 

numbers recaptured per 
unit effort (y, ) 

Log ( yi /T) 
Weighting factor 



1-12 

1 

139 
-1.974 

0.139 



13-24 25-36 37-48 

234 

91 52 40 

-2.397 -2.957 -3.219 



0.091 



0.052 



0.040 



T - 1000 

Weighted regression of log (y,/T) on i gives 
log (y t /T) - -1.539 - 0.438i 
Hence Zp - 0.438 
Hence Z - 0.438 
S - exp -Z 
- 0.645 
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Tbble 7 

Estimation of survival rate from a single tagging experiment 
given exact recovery times for each fish recaptured 



Numbe r o f days absence for i n JL 2 '*5tL r<-> oa^ tore ri. 

1, 5. 7, 9, 10, 12, 16, iS, ;/-J. : r , v. s 30, 50 f 7 T 7 c ?r ^r *\ , H , 162, I9i> 25 
Mem duration of fl^eedoiii t ..-.. ;->, O^VT. () I <' years 
Oulland (1955): 

Z( annual) - 1/t - i^ * r ;, 

u ,' * 

Paulik (1963): 

t 0.17 

u - 250/365 m 0.68 year 
t/u - 0.25 

Henoe Zu (which is equivalent to o'x in Table 1 of Deemer and Votaw, 1955) 

- 3-55 
/. Z 3*55/0.68 . 5,2? 



Table 8 

Estimation of survival rate 
Jolly Se"ber method 



Sampling times 


Numbers ? - 
tagged (T,,) * * 


456 


r i 


1 


1 000 


60 32 


21 13 3 


129 


2 


2 060 


124 


82 49 12 


267 


3 


1 656 




149 89 22 


260 


4 


1 252 




188 46 


234 


5 


2 339 




193 


193 




m. 
i 


60 156 


252 339 276 


1 083 


Values 


of q 








1 


2 


3 4 


5 







69 


180 188 


83 




Values 


of S 








S 


1/2. 


129 [(2060) (69) + 
0.59 


(267X60)]/(267)(1 000)(0 


+ 129) 


S 


2/3- 


267 [(1 656)(180) 
0.50 


+ (260X156)]/(260)(2 060)(69 + 267) 


S 


3/4 - 


260 [(1 252)(188) 
0.45 


+ (234)(252)y(234)(l 656) 


( 180 + 260) 


S 


V5- 


234 [(2 339) (83) H 
0.60 


H (I93)(339)y(l93)(1 252) ( 


188 + 234) 
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Table 9 

Estimation of survival rat - Robson method 
Samples for tagging and recaptures taken during discrete intervals 



Numbers 


Sampling 


times 






2 3 


4 


r. 


1 1 000 
2 2 OOO 
3 1 500 


60 36 
120 


27 
90 

135 


123 
210 

135 


m 60 156 


252 


Values of q: 




1 2 3 
63 117 




Values of Si 




Sy 2 - (2 000)(63)(123)/(1 
- 0.60 


000)(123)(210) 


S 2 ^ - (1 500)(l17)(210)/(2 
- 0.50 


000)(63 + 210)(135) 
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Table 10 

Estimation of survival rate - Rob son method 
Recaptures made continuously between taggings 



Recapture periods 


,v Numbers _ 2 
V ; tagged (0^) ^ J 


3-4 4-5 5-6 


r i 


1 1 OOO 100 72 


60 34 12 


278 


2 2 000 240 


200 112 40 


592 


3 1 500 


300 169 61 


530 


4 1 000 


250 90 


340 


5 


300 


300 


m 100 312 


560 565 503 


2 040 


Values of q 






1234 


5 




178 458 428 


203 




Values of S 






S 1/2 . (2 000)(178)(278)/[1 
- 0.60 


000(0 + 278)592] 




S 2/3 - (1 500)(458)(592)/[2 
0.50 


000(178 + 592)530] 




S 3/4 - (1 000)(428)(530)/[1 
- 0.45 


500(458 + 530)340] 




S 4/5 - (2 000)(203)(340)/[1 

m 0.60 


000(428 + 340)300] 
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Table 11 
Data from Table 9 for application of Bailey's triple catch equation 



Tagging 


time 


Numbers 
tagged 


Recoveries at 
recapture times 

2 3 


1 




1 000 


60 36 


2 




2 000 


120 


T 1 - 


1 000 
? 000 


3 1/2 " 
3 1/2 * 


? 000(36) , 




i 000(120) " ' b 



Table 12 

Data from Table 8 for application of Leslie and Chatty's 
triple cntch ecjuatiori 



Tagging time 


Numbers 
trigged 


Recoveries at 
recapture times 

2 3 


1 


1 000 


60 32 


2 


2 000 
i 


124 


60 156 


T 1 - 1 OOO 
T 2 .2 000 


S 1/2 - f(2 000) (32) + (156)(60)]/(1 000)(124) 
- 0.59 
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Table 13 
Estimation of P from dnta in Tables 4 and 6 



From Table 4: 



Intercept I - 5347 
T - 1 000 



P m 0.430[exp(5.347 - 0.430)1/1 000(l-exp( -0.430)) 
.0.17 



Prom Table 6: 



Intercept I - -U539 

T - 1 000 



F - 0.438[exp(--1. 539-0.433) V(l-e:rp( -0.438)) 
,0.17 



Table 14 

Estimation of F from data in TVible 5 
Gulland method 



Intercept I 



> 132 



T - 1 000 

F - (exp 5.132)/1 000 
- 0.17 



Table 15 
Etimation of P from data in Table 7 



T m 1 000 

m - 21 

"t - 0. 17 years 

Gulland 1 s method: 

F - 21/1 000 (0.17) 
- 0.12 

X - Z - P - 5-9 - 0.12 - 5.78 

Chapman's method: 

P - 20/1 000 (0.17) 
. 0.12 
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Table 16 
Ai Estimation of the rate of capture (p) using tho data in Tnble 8 





T i 


Baaio 


data 


m. 
i 


Estimated 


values 


i 


r. 


h 


1 


1 000 





129 


_ 


~. 


M 


2 


2 060 


69 


267 


60 


592 


0.10 


3 


1 656 


180 


260 


156 


1 302 


0.12 


4 


1 25^ 


188 


234 


252 


1 258 


0.20 


5 


2 339 


83 


193 


339 


1 345 


0.25 



B: Estimation of the rato of capture (p) usin the data in Tnble 9 







Basic data 




Estimated values 




T, 




m 


I' 




i 


1 L 


i 


^ x 


1 


1 OOO 


123 








2 


2 000 


63 210 


60 


600 o. 10 


3 


1 500 


117 135 


156 


1 300 0.12 



V - h - 
Pi - "i/I 



Table 17 
Estimation of F using the data in IVble 10 



i 


T i 


* 


r i 


m i 


s i 


z i 


H i 


P i 


1 


1 000 





278 


100 


0.60 


0.51 


1 000 


0.13 


2 


2 000 


178 


592 


312 


0.50 


0.69 


2 601 


0.17 


3 


1 500 


458 


530 


560 


0.45 


0.80 


2 796 


0.29 


4 


1 000 


428 


340 


565 


0.60 


0.51 


2 559 


0.28 



F . 



1 - 3) 



34 



FUS/T153 



Table 16 
Formulae for estimating population size (IT) by the Peterson method 





T^pe of 


Bi 


Btimates of 


Reference 


Sampling 


Population size (N) 


Variance of (H) 










ATUilnv 1QS1 


TH **+ 


K-^2 


ITAT IT ^ ^ n ( n "" m ^ 


BHa4 1 Air 1QC1 1 


T\4 - r>A r + 


N " IT 

N T(n 4- 1) 


in 
O^Cn + 1)(n - m) 


JtKtXXejTf \Jj I 

CPVkAYMnMVi 4QCM 


TVi fw^/>4* 


11 " m + 1 
(T 4- 1)(n + 1) 


(m + 1) (m + 2) 
var H - B 2 /^ + 2( S ) 2 + 6( 2 )3\ 


Sohaefer, 1951 

"H fTt4"l*tr 1OR^ 




11 - m 4- 1 - 1 
n(T 4- 1) 


^"^ F ' fl InT * A nT ; + b ^ nT ; J 
v (T - m + 1)(E + 1)(K - T) 


u* iJBalieyt liOv 
(Chapman, 1952) 


Inverse 


m 


v ~ ff - m(T+ 2) 



Estimation of population size by the Jolly and Sober method 



i I, ^ n, H, 


1 





^p 


MB * 





2 


592 


60 


16 000 


158 ooo 


3 


1 302 


156 


10 000 


84 ooo 


4 


1 258 


252 


7 ooo 


35 ooo 


5 


1 345 


339 


5 500 


22 000 


Values 


of I 4 


and m, from 


Table 16A 






i 


i 






Values 


of n. 


supplied for purposes of the 


example 
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Table 20 
Example of the estimation of growth parameters from tagging dat* 



a 


Lt 


L t + a 


y 


X 


b 1 






Time of 
liberty 
(years) 


Length at 
release 
(cm) 


Length at 
recapture 
(cm) 


L+ ^ - Lt 
t 4- a 


L. _,_ -I- Lt 
t 4- a 


i* a 


\ 


b 1 


a 


2 


tanh b 


y 'tanh b 1 


0.5 


9-1 


16.9 


15.6 


13 


0.05 


1.0008 


15.6 


0.8 


24.2 


33.8 


12*0 


29 


0.08 


1 . 002 1 


12.0 


1.0 


41-5 


50-5 


9.0 


46 


0.10 


1.0033 


9.0 


0.4 


61.9 


64.1 


5-5 


63 


0.04 


1.0005 


5-5 


1.2 


74.2 


77.8 


3.0 


76 


0.12 


1.0048 


3.0 



Table 21 
Example of the calculation of movement parameters from tagging data 



p (km) 


e 4 


todays) 


sin e i 


cos 8. 


r.sin 9. 


r.cos 8, 


X 


5 


58 


18 


4-.84B 


*-530 


f 4.24 


4- 2.65 


1.4 


50 


150 


60 


f.500 


-.866 


f25.00 


-43.30 


41.7 


20 


230 


30 


-.766 


-.643 


-15-32 


-12.86 


13.3 


30 


335 


45 


-.423 


+.906 


-12.69 


4-27.18 


20.0 


Totals 




153 






+ 1.23 


-26.33 


76.4 






R - V 1.23 2 + (-26. 33) 2 - 26.36 


V - 26.36/153 - 0.17 km/day 


2 1 <-7* A 694.3, , A . 2 7 , 
a - -T (76.4 ifj"! " 18 km /^y 
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Table 22 

Showing some possible combinations of values of c and v 

P 

such that ev 4*0 km /day 



(km) c 


v (kn/day) 


4 


1 


1 


4 


0.1 


40 


0.01 


400 



PIRS/T1S3 



1 



Showing a way of classify ing tag returns 
Sampling Times 



2 3 


456 


7 


Total rcapturB 


"12 m 13 


m l4 m l5 m 16 


"17 


r 1 


"23 


n 2A m 25 26 


m 27 


T 2 




m M m 35 m 36 


m 37 


r 3 




m 45 m 46 


"47 


r 4 




"56 


m 57 


r 5 






"67 


r 6 



Total recaptureo? 



m. - number liberated on the ith sampling occasion and recaptured on 
3 the jth sampling occasion; 

m, - total number recaptured on the jth sampling occasion; 
J 



r. - number recap-tared out of the fish liberated on the ith sampling 
occasion* 



-3L 
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Fig. 2 ESTIMATION OF SURVIVAL RATE 
Method of Jolly and Seber 




For the 3rd occasion 

m 3 )= r 3/ /T 3 
.M 3 = m 3 + q 3 T 3 /r 3 



m 3 



>r 3 



For the /.th occasion 



m. 



The proportion surviving from the end of the 3rd to the end of the 
occasion is given by: 



In general 



r t T 3 (q 3 *r 3 ) 
r t r 1 "litJL 



+ 1 T,(q, 
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Fig. 3 ESTIMATION OF SURVIVAL RATE 
Method of Robson 



ccasion 



= T 



J 



m: 



For the 4th occasion 



The proportion surviving from the end of the 3rd to the end of the 4th 
occasion is given by: 



In general 



_ Tt + 1 qt + 1 r> 
T t r t ^ (q t ^r t ) 
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Fig. I* ERRORS OF TYPES A, B & C 



in number 



regression in the absence 
of error 

regression with error 



TYPE 'A' ERRORS 

influence the intercept 

but not the slope of 

the regression 



time 



In number 






\ 



time 



TYPES 'B'&'C' ERRORS 

influence the slope but 

not the intercept of 

the regression 



FIR3/T153 



Fig. 5 SHOWING THE RELATIONSHIP BETWEEN THE 
POSITIONS OF LIBERATION (A) AND 
RECAPTURE (B) OF A FISH WITH REFERENCE 
TO THE CO ORDINATE SYSTEM NESW 



N 

A 



w- 



B 



S 



Fig. 6 



B 




Fig. 7 DETERMINATION OF THE OVERALL GROUP 
DISPLACEMENT, R 



W- 



N 



y = I r cos 6 




x = Ir sin 6 



R 



-11185 



